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ABSTRACT 

We investigate the X-ray absorption structure of oxygen in the interstellar 
medium by analyzing XMM-Newton observations of the low-mass X-ray binary 
Sco X-1. We use simple models based on the O I atomic photoabsorption cross 
section from different sources to fit the data and evaluate the impact of the 
atomic data on the interpretation of the observations. We show that relatively 
small differences in the atomic calculations can yield spurious results, and that 
the most complete and accurate set of atomic cross sections successfully reproduce 
the observed data in the 21.0—24.5 A wavelength region of the spectrum. Our 
fits indicate that the absorption is mainly due to neutral gas with an ionization 
parameter of ^ = 10~^ erg cm s~^ and an oxygen column density of Nq ~ 
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(8—10) X 10^^ cm~^. Our models are able to reproduce both the K edge and the 
Ka absorption line from O I which are the two main features in this region. We 
find no conclusive evidence for absorption by other than atomic oxygen. 



1. Introduction 

X-ray spectroscopy provides a powerful tool for understanding the physical and chemical 
properties of the diffuse interstellar medium (ISM). The X-ray band covers the emission 
and absorption spectra produced by inner-shell transitions of the most abundant ions from 
carbon to iron. The interaction of X-ray photons from bright point sources (e.g. galactic 
X-ray binaries) with the ISM gives rise to absorption lines and edges in the spectrum. The 
energy position and the shape of these features depend on whether the absorption is due to 
free atoms or molecules, and on whether these atoms or molecules are in the gas or in the 
sohd phase. 

Neutral oxygen is a major constituent of the ISM which makes it one of the most im- 
portant elements in astronomical observations. Precise knowledge of the neutral oxygen 
atomic parameters is needed for the correct modeling of the observed spectra. Calcula- 
tions of the photoabsorptio n cross section of the ground state of O I were carried out by 
McLaughlin &: Kirbyl ( 119981 ) (hereafter MC K98) using th e .R-m atrix method, giving a detailed 
comparison with the experimental results of lStolte et al.l ( 119971 ). Although they claimed over- 
all agreement, there are significant discrepancies in the positions of the inner-shell excited 
resonances and in the nea r-threshold resonance pro f iles. T his problem was overcome in the 
LS'-coupling calculation of iGorczyca fc McLaughlinI ( l2000l ) (hereafter GMCOO) by taking into 
account core relaxation effects and the smearing of the K edge due to Auger damping. A 



more complete matrix calculation was carried out in intermediate coupling by iGarcfa et al. 



( 120051 ) (hereafter gar05) for all the ions in the oxygen isonuclear sequence. There is now very 
good agreenaent b etween the GAR05 calculations and both the experimental cross section 
( jStolte et al.lll997l ) and the GMcOO results. 



The oxygen inner-shell features in the X-ray spectrum of galactic sources have been 
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2OO2I : iTakei et al.ll2002l : I Juett et al.l 12004 : iTurner et al.l 12004 : 



20061 ). Howe ver, studies in the IR and U V have shown that 



oxygen can also be found in solid particles ( iDraind 120031 : IWhittetl 120031 ). It has been argued 
that oscillatory modulations near the K edge, usually referred to as the X-ray absorption fine 
structure (XAFS), could be detected. These are co ndensed matter modulations of the atomic 
cross section due to the presence of solid particles (|Lee &: Ravelll2005t iLee et al.ll2009l ). Stud- 
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ies of the soft X-rays from galactic sources have reported possible detections of molecules 
that could be linked to XAFS signatures in the edges of sev eral elements such as Ne, Si, and 



Mg JPaerels et al.lboOll: Ihee et all [20021 lUeda et al.lboosi ): in the L edge of Fe Jlee et al. 



20011: 



Kaastra et al 



200 



and in particular, in th e oxygen K edge (jPe Villiers et al. 



2003 



Costantini et al.l l2005l : Ide Vries fc Cost ant inil l2009l : IPinto et al.l I2OIOI ) . Although these sig 



natures could be important in the K edges of molecules involving higher Z elements, namely 
Mg, Si, and Fe, oxygen is ~ 10—20 times more abundant than any of these, thus poten- 
tially providing enough signal-to-noise to detect the XAFS signatures. Theoretical model s 
of X AFS in the astrophy sical cont ext have be en developed by lWod ( 119951 ) : IWoo et al.l ( 119971 ): 
and lForrey et al.l ( 11998! ) . See also iLeel (120 111 ) and references therein for details on the theory 
of XAFS. 

XMM-Newton observations of the X-ray source Scorpius X-1 (Sco X-1) reveal strong 
absorptio n in the wavelength reg ion corresponding to neutral oxygen. Located at a ~ 2.8 kpc 



distance (IBradshaw et al.l Il999l ) and with a flux of F ~ 3.4 x 10^ erg cm ^ s ^ (in the 



2—10 keV energy band), it is the brightest X-ray source in the sky other than the Sun and 
the diffuse X-ray background radiation. Its high X-ray flux provides very good statistics 
in relatively short exposure times, giving the opportunity to study signatures of oxygen 
absorption in the IS M with great detail. For several galactic sources, including Sco X-1, 
de Vries et al.l ( 12003! ) analyzed high-resolution X-ray spectra taken with the reflection grating 
spectrometer (RGS) in the XMM-Newton satellite. By comparing low and high extinction 
sources, they w ere able to separate the ISM and the instrumental components of the O I K 
edge; moreover, Ide Vries &: Costantinil (120091 ) searched for XAFS signatures in the spectrum 
of Sco X-1. The XAFS signature is derived from the differences between the observed flux 
and that predicted theoretically. However, the model used by these authors is based on the 
atomic oxygen absorption cross section calculated by mck98. 

In this Letter we show the importance of the accuracy of the atomic data used in the 
modeling of the detailed features of the oxygen absorption in the ISM. We demonstrate 
that small variations in the K-edge structure derived from different atomic calculations yield 
spurious results when applied to astronomical observations. In Section [2l we describe the 
observational data used in our study while the theoretical models are delineated in Section 121 
Results derived from model fits of the observed data are presented in Section HJ and finally, 
the main conclusions are summarized in Section [51 
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2. Observation and data reduction 



For the purpose of this work, we adopt the XMM-Newto n spectrum of Sco X-1 taken 

in orbit 0592 (Obs ID = 0152890101) with the RGS i nstrument Jden Herder et allboOlh. The 

genera l observational strategy is described in detail in lde Vries et ahl (|2003[ ) and lde Vries fc Costantini 
( 120091 ) and will not be given here. Taking into account the high X-ray flux of Sco X-1, stan- 
dard spectroscopy mode would lead to a very high pile-up level. This was avoided by choosing 
a faster readout mode which is able to read each of the nine RGS CCDs in separate expo- 
sures. Because of a failure in the RGS2, we make use of only the RGSl data in what follows. 
During this observation there were 24 separate contiguous exposures, but only the fifth and 
the eighteenth contained calibrated and reliable data in the 21—25 A spectral range; thus, 
we only use these exposures in our spectral analysis. In the notation of the SAS rgsproc 
task, these two exposures are labeled S005 and SO 18. We follow the standard procedure for 
the reduction and extraction of the RGS spectrum using SAS version 10 with the latest 
calibration files (CCFs). We finally present the spectrum with the default binning of 0.05 A 
employing full spectrometer resolution. 



3. Spectroscopic models 
3.1. Atomic data 

Figure [1] shows the photoabsorption cross section of the ground state of neutral oxygen 
in the 22.5—24.0 A wavelength region from the calculations by mck98, GMcOO, and GAR05. 
All these curves have been convolved with a Gaussian profile of 0.182 eV full width at half 
maximum (FWHM) in order to match the resolution of those presented in GMcOO. This 
spectral region covers both the K edge around 22.5 A and the Ka Is— 2p absorption line 
near 23.5 A. Besides the differences in the Ka position, the GAR05 and GMcOO calculations 
agree very well. Nevertheless, there are significant discrepancies with MGK98, particularly 
in the shape of the K edge and in the energy separation between the K edge and the Ka 
resonance. Since the Ka absorption feature is prominent and well-resolved and its energy is 
experimentally fixed, it is used as a reference in most spectral fits. Thus, uncertainties in 
the energy and shape of the cross section around the inner-shell edges are translated in the 
form of spurious residuals in the fits. 

These discrepancies have been formally addressed by GMcOO as the result to two con- 
spicuous effects: orbital relaxation due to the K shell vacancy which affects the atomic 
structure, and consequently, the resonance and edge positions; and resonance decay through 
an infinite number of transitions wherein the resonance width is usually dominated by the 
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spectator- Auger channel. It is evident from the comparison in Figure [T] that the neglect of 
these two effects underestimates the absorption cross section near the K edge by 20—50%. 



3.2. Model A 



In order to resolve the features in the Sco X-1 observation, we have employed the 
photoionization code xstar. Several calculations were performed covering a range of pa- 
rameters, the most important being the hydrogen column density, N-^, and the ionization 
parameter, ^ = L/nB?, where L is the lu minosity of the source, R its distance, and n is 
the density of the gas (ITarter et al.l 119691 ) . We have constructed a grid of xstar models 



spanning hydrogen column densities of Ah = 10 —10 cm~ and ionization parameters of 
^ = 10~^— 10 erg cm~^ s~^ with the gas density fixed at n = 10^ cm~^. The spectral region 
of interest, 21.0—24.5 A, is relatively small and, in practice, is dominated by oxygen species. 
Thus, the xstar models include hydrogen, helium, and oxygen ions in the ionization bal- 
ance calculation, assumin g an oxygen abundance relative to hydrogen of Aq = 6.8 x 10~^ 
(iGrevesse &: Sauvall Il998l ). We will refer to this model as Model A. The xstar models 
incorporate the GAR05 cross sections for all the oxygen charge states. 



3.3. Models B, C and D 

If a simple description of X-ray absorption by a cold medium assumes that all the 
spectral features are solely due to neutral oxygen photoabsorption, the observed flux can 
then be approximated as 

F{E)=Foexp[-No,ao,{E)] (1) 

where Fq is a normalization factor, A"oi is the oxygen column density, and croi(-E') the pho- 
toabsorption cross section of O I. This is a convenient way to evaluate the relevance of the 
atomic data used to fit the observation. Using this assumption, we have produced three 
additional models using Equation ([1]), each with a different cross section: Model B includes 
that computed by MCK98; Model C by gar05; and Model D by GMCOO. Note that the 
XSTAR model (Model A) and Model C are equivalent in the sense that they use the same 
atomic data although the former also includes the background due to H and He. These 
models are summarized in Table [H 
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4. Results 

The models described in Section [3] are used to fit the absorption features observed in the 
X-ray spectrum of Sco X-1. In this respect and to determine the corresponding statistics, 
we use the X-ray spectral package XSPEC vl2.3.0, and all the fits are carried out in the 
21.0—24.5 A spectral region. Figure |2] shows the fit using Model A, our main model, since it 
is the result of a self-consistent photoionization calculation and incorporates the most recent 
atomic cross section for O I by GAR05. In the upper panel, the black and gray data points 
belong to exposures S005 and S018, respectively. The best fit using Model A is shown with 
solid fines, red and blue respectively corresponding to the fit applied to each exposure. The 
spectrum displays a very strong atomic K edge which covers the 22.8-23.3 A wavelength 
range. The Ka absorption resonance at ~ 23.5 A is also one of the strongest features while 
the K/3 is much weaker but is still detectable at ~ 22.9 A. In the lower panel, we show the 
residuals with respect to the model in a units. The black and gray points correspond to 
exposures S005 and S018, respectively. Model A fits the K edge and the Ka absorption line 
successfully in the two exposures of Sco X-1. Statistics for the combined fit (i.e., for the 
two exposures) show a reduced chi-square of x^/dof = 2.75620 (where dof is the number of 
degrees of freedom). The best-fit hydrogen column density is A^h = (1.33 ±0.02) x 10^^ cm~^ 
which corresponds to an oxygen column density of Nq = (9.04 ± 0.12) x 10^^ cm~^. 

In order to get a better grasp of the atomic data effects on the description of the observed 
spectra, we now fit the measurements with the simple models based on the raw photoabsorp- 
tion cross section of neutral oxygen (see Equation [1]). In Figure [3] we depict the fits using 
Models B and C; in the upper panel, the black/gray data points are the measurements while 
the red/blue curves are the models corresponding to exposures S005/S018, respectively. The 
best fits using Models B and C are respectively shown with dashed and solid lines. The mid- 
dle panel gives the residuals in a units with respect to Model B and the lower panel those 
with respect to Model C. It is important to note that, if these two models were equivalent, 
the dashed and solid lines with the same color would be close to each other. However, there 
is a clear discrepancy between the dashed lines (Model B) and the solid lines (Model C) in 
the region near the oxygen K edge (~ 22.5—23 A). These differences may be appreciated 
in the residuals of the Model B fit shown in the middle panel. Furthermore, the residuals 
also indicate that Model B cannot completely fit the intensity of the Ka absorption line at 
~ 23.5 A; it predicts an oxygen column density of A'o = (9.25 ± 0.07) x 10^^ cm~^; and the 
fit statistics give a reduced chi-square of x^/dof = 6.02484 which is poorer than Model A. 
As expected, the fit using the raw cross section by GAR05 (Model C) is equivalent to that 
using Model A giving a reduced chi-square of x^/dof = 2.44650. Model C predicts an oxygen 
column density of A^o = (7.94 ±0.22) x 10^^ cm^^ somewhat smaller than Models A and B. 
The differences between Models A and C may be due to the numerical interpolation used in 
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the storage and retrieval of the cross sections by the xstar package. 



Model B is the same as that used by Ide Vries &: Costantinil (120091) to fit the oxygen 



absorption in the spectrum of Sco X-1. These authors claimed the detection of XAFS 
signatures in the spectra based on the relative changes in the observed flux with respect to 
the smooth flux predicted by the model for energies above the K edge. They also argued 
that the apparent shift of the observed edge with respect to the atomic model could be due 
to the fact that some fraction of the oxygen in the ISM is bound in solids. However, the 
analysis presented here shows that the large residuals found in the Model B fit are artifacts 
brought about by the mck98 cross section. 

Finally, in Figure|l]we show a similar comparison using the GAR05 and the GMCOO cross 
sections (Models C and D). As before, the upper panel shows the Sco X-1 spectrum and the 
best-fit models. The black/gray data points are the observations while the red/blue curves 
are the models corresponding to exposures S005/S018, respectively. Note that Figure H] 
covers the 22.5—24.0 A wavelength range to enhance the K-edge region. The best fits using 
Models C and D are respectively shown with solid and dashed lines. The middle and lower 
panels give the residuals in a units for Models D and C. The two models are equivalent giving 
similar fits (x^/dof = 2.76445 for Model D). This is consistent with the accord between the 
GMcOO and GAR05 cross sections. The oxygen column density derived from this model, 
Nq = (10.49 ± 0.06) X 10^'' cm~^, is larger than those estimated by all the previous models, 
but consistent within the uncertainties in the parameters predicted by Models A and C. The 
discrepancies in the predicted oxygen columns in Models C and D can be explained by those 
in their respective cross sections. In this small spectral region, the oxygen column depends 
almost entirely on the depth of the K edge. By comparing the values of the cross section at 
22.8 and 23.0 A (i.e. before and after the K edge) in Figured! it may be noticed that the K 
edge in GMcOO is weaker than in gar05. If a feature is weaker in a model, it will require 
a larger column density to reproduce the observation and, therefore, explains the different 
column densities obtained. All the proposed models and the parameters derived from their 
corresponding fits are summarized in Table [TJ 

In all the fits presented here, we see residuals of around ±2cr distributed homogeneously 
along the entire spectral range considered in our analysis, consistent with the reduced chi- 
square close to 2 obtained in our best fit. This suggests that most of the fit errors are 
systematic. However, we notice significant residuals at wavelengths shorter than the Ka line. 
Due to a gap in the 23.33—23.44 A wavelength range, the detection of spectral features is not 
possible; nonetheless, both the atomic Ka resonance due to ionized oxygen and the Is n* 
resonance due to m olecular oxygen oc cur at 23.35 A. The former alternative is consistent 



with the analysis by lJuett et al.l (120041 ) where they detected O ll absorption in the spectra 
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of several X-ray binaries. On the other hand, our fit using Model A predicts that absorption 
occurs in a mostly neutral gas with an ionization parameter of ^ = 10~^ erg cm s~^ which 
would also suggest molecular absorption. Therefore, given the lack of data in this particular 
spectral range, we cannot rule out the presence of molecules in the observation in hand. 

5. Conclusions 

In the present report we have shown the relevance of accurate atomic data in the mod- 
eling of the X-ray spectra from cosmic sources. In particular, wc have studied the Sco X-1 
spectrum produced by the RGSl instrument on board of the XMM Newton satellite cover- 
ing the 21.0—24.5 A wavelength region. Absorption occurs when the X-rays interact with 
the cold gas of the ISM, the main spectral features in this region being the absorption K 
edge and Ka line from neutral oxygen. We found a good fit using a self-consistent pho- 
toionization model which includes the most recent atomic data for the oxygen isonuclear 
sequence by GAR05. Our fits indicate that the absorbing gas has an ionization parameter of 
^ — 10~^ erg cm s~"^ and a hydrogen column density of A^o ~ (8—10) x 10^^ cm~^. 

Simple models based on the raw atomic photoabsorption cross sections of O I from 
three different calculations were used to evaluate data sensitivity. We show that models 
based on the mck98 atomic cross sections are unable to reproduce the K-edge structure in 
detail, while those based on GMCOO and GAR05 yield more accurate fits of the main spectral 
features. The fits using the most up to date models do not show evidence for absorption by 
anything other than atomic oxygen. 

The analysis presented here indicates that the atomic data uncertainties in combination 
with the limited resolution of the grating spectrum make detection of molecular or solid 
material challenging. Although oxygen is expected to be found in molecular form or locked 
into solids in the ISM, the use of accurate atomic calculations to correctly account for the 
atomic oxygen contribution is crucial when searching for XAFS or similar features in the 
X-ray spectra of astronomical sources. 

We thank Tom Gorczyca for providing the MCK98 and GMCOO calculations. This work 
was supported by a grant from the NASA astrophysics theory program 05-ATP05-18. This 
research has made extensive use of the NASA Astrophysics Data System. 
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Fig. 1. — Comparison of the theoretical pho toabsorption cross sec t ions o f neutral oxy- 
gen in the 22.5—24.0 A wave length region: McLaughlin fc Kirbv ( 19981 ) (blue curve); 
Gorczyca fc McLaughlinI (120001 ) (red curve); and iGarcia et al.l ( l2005l ) (black curve). This 
spectral region covers both the absorption K edge and the Ka absorption line (Is— 2p) from 
O I. All the curves have been convoluted with a 0.182 eV FWHM Gaussian profile. 
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Fig. 2. — The Sco X-1 spectrum produced by RGSl aboard XMM-Newton covering the 
21.0—24.5 A wavelength region. In the upper panel, the black and gray data points corre- 
spond to exposures S005 and SO 18, respectively, while the best fits to each exposure using 
the XSTAR photoionization model (Model A) are respectively depicted in red and blue. The 
lower panel shows the residuals in a units in the same range with respect to the model. 
Black and gray points correspond to exposures S005 and SO 18, respectively. 
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Fig. 3. — Spectral fit comparison using two different models. Upper panel: Sco X-1 spectrum 
covering the 21.0—24.5 A wavelength region. The black/gray data points are the observed 
data while the red/blue curves are the models corresponding to exposures S005/S018, re- 
spectively. The best fits using Models B and C are shown with the dashed and solid lines, 
respectively. Middle panel: residuals in a units with respect to Model B. Lower panel: 
residuals in a units with respect to Model C. 
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Fig. 4. — Spectral fit comparison using two different models. Upper panel: Sco X-1 spectrum 
in the 22.5—24.0 A wavelength region. The black/gray data points are the observed data 
while the red/blue curves are the models to exposures S005/S018, respectively. The best fits 
using Models C and D are shown with solid and dashed lines, respectively. Middle panel: 
residuals in a units with respect to Model D. Lower panel: residuals in a units with respect 
to Model C. 
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Table 1. List of Models 



Model 


Description 


Atomic Data 








Notes 








(10" cm-2) 








A 


Full XSTAR model 


AKAOI^ +GAR05'' 


9.04 ±0.12 


2.75620 


No 


significant residuals 


B 


Atomic cross section'^ 


mck98'^ 


9.25 ±0.07 


6.02484 


Lar 


ge residuals near K edge 


C 


Atomic cross section'^ 


GAR05'' 


7.94 ±0.22 


2.44650 


No 


significant residuals 


D 


Atomic cross section'^ 


GMCOO'* 


10.49 ±0.06 


2.76445 


No 


significant residuals 



HBautista fc KallmanI l|200ll ') 
iGarci'a et Zl | |2005| '| 
' iMcLaughlin fc Kirbvl lll998l ') 
' ^Gorczvca fc McLaughlinI hoodl) 
^See Equation |(T]| 



